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Onc-dimensional problem for quantum harmonic oscillator with "regular+random" 
frequency subjected to the external "regular+random" force is considered. Aver- 
aged transition probabilities are found. 



1 Introduction 

There are a lot of works published recently and devoted to quantum chaos, 
i.e. they consider quantum systems analogous to classical systems that demon- 
strate the features of chaotic behaviour. Exploration is developed at different 
directions: investigation of energy levels distribution; derivation and calcu- 
lation of objects (analogous to classical Lyapunov exponents or KS-entropy) 
which can witness for chaos in quantum system and others. 

We present a new approach for quantum description of above mentioned 
systems that may be used for the investigation of a problem of multichannel 
scattering. It is known that scattering with rearrangement (typical for chemical 
reactions) goes through resonance complex formation which causes chaotic 
behaviour. It is impossible to predict the way to be followed by a system after 
leaving the complex because of the small parameters change. As it was shown 
in our previous works, Bu the scattering process with rearrangement may be 
described in framework of randomly walking harmonic oscillator model. In this 
paper we consider_one-dimensional case. More simple case has been already 
introduced earlier 

2 Description of the problem. 

We consider the wave function ^ stc (t, x), describing the state of the system as 
a random process with the time evolution determined by equation 

id t * stc = HV stc , (1) 

where ID Hamiltonian H is quadratic in space variable 

* = -\l^ + \ n2(t)x2 - nt)x ' (2) 



while functions f2 2 (t) and F(t) are random functions of time variable. We 
suppose that by definition 



n 2 (<) = ng(t) + v^/iWefi - ti), 



(3) 



where and -Fo(^) are deterministic functions while /2(f) are inde- 

pendent zero mean gaussian random processes with two-point correlations of 
5— form: 

< f i (t)f j (t')>=S ij S(t-t'), i,j = 1,2. (4) 



Constants i 



1, 2 control the power of forces fi(t), t — 1,2, while functions 



Pi(t), i = 1,2 arc nonnegative Px,P2 > 0, which leads to the suggestion of the 
following asymptotic behaviour 



fio(t) 



— > n/,» v F (t) ■ 

f — >±oc lout; r— >±oo t — >-±oo 



0, Pi(t) 



0, i = 1,2. 



(5) 



which guarantees in the limit t — > — oo existence of stationary states <p l "(t,x) 
<j>™(t,x) = e- i ( n+1 /V n *» t <f>™(x), 
( I — — \ 1/2 

= I airV^ 1 e-^^/^^yn-^), (6) 



where H n (x) are Hermitian polynomials. At the limit t — > +oo there also exist 
stationary states </>™'(i, x), which may be got from (^) by simple replacement 
of Qi„ by flout- Moments t\ and £2 of switching on the noise are chosen to 
be finite in order to provide the correctness of subsequent constructions. Now 
we formulate our task: to get the averaged transition probabilities W nm from 
initial stationary states 4> l "(t,x) to final ones x) which come from the 

evolution described by (|l|)-(0). 



3 Formal expressions for wave functional and transition probabili- 
ties. 

The main results of this paper are found on the basis of formal solution of 
(|l])-(||), which may be constructed for arbitrary Q 2 (t) and F(t) as a functional 
of solutions of classical equations of motion. Namely, as was shown inta, the 
following representation for solution of (^)-(^) may be written that turns at 
the limit t — > — oo to the stationary state <j>™(t, x) defined in (||): 



-.exp< i 



rj(x-rj) + — (x-rj) 2 + a 



4%b 



'), n = l,2,. 

(7) 
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where 77(f) is a solution of the classical equation of motion for oscillator with 
frequency f2(f) subjected to external force F(t) : 

ij + n 2 (t) v = F(t), 7?(-oo) = tK-oo) = 0, (8) 

tr(f) is an action functional corresponding to this solution 



a(t) = 



If, 2 - -fiV + Fr, 



dt', (9) 



while r(f) and r(f) can be expressed in terms of solution of homogeneous 
equation corresponding to (|J) 

e + n 2 (f)£ = 0, £(t) ~ e iQ ^ (10) 

t—> — OO 

in the following way: f(t) = r(f)e 7 W, r(t) = |£(f)|, r = j(t)/Q in . Func- 
tions ^>^l(t,x) are in fact functionals for the realisation of random processes 
fi(t) and /2(f)- Therefore it is natural to call them "wave functionals". 

We are interested in construction of transition probabilities from the states 
^s^l (f, x) to stationary states </>™ *(f, x) in the limit f — > +00 averaged over the 
random processes /i(f) and /2(f) realisations. Let's designate them as W nm . 
Having defined 



V%l(t,x) = ^c nm (t\fi,f 2 )(t>™ t (t,x), W nm = lim <|c„ m | 2 ), (11) 

m=0 



we denoted the averaging over /1 and /2 by the symboL,< .. > . An expression 
for generating function of coefficients c nm is known (seeB). Having at hand for- 
mal expressions for objects of interest we should turn the averaging procedure 
to the form convenient for analytical or numerical treatment. 

4 Equation for the distribution function. 

Functions £(f) and 77(f) defined by equations (JsJ) , ( |l0| ) are random processes. 
It is more convenient to treat with the other random processes z% = 77, 22 = 
77, Z3 — -Re (£/£), Z4 = Im(£/£) which in total we denote briefly as z. One 
can find that the Focker-Plank equation for conditional distribution function 

P(z,t\C,t>)= (S(z(t)-z)) 



, p looks as follows: 

z ( t > )—i 



dP J^d(KiP) , 2N a 2 P d 2 P n d 2 P - 

-m=-l i SzT + {e2P2 +eiPlZl) ^zJ +eiPl ~dzJ + 2eiPlZl dSz~ 3 = LP > 

(12) 
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where K\ = z 2 , K 2 = F — OqZi, K 3 = z\ — z\ — Qg, K4, = — 2z 3 z 4 . It 
is supplemented with the initial condition P\t=t > = 5(z — Q and with the re- 
quirement that J Pdz is finite. Random vector £ of initial values of trajectories 
is described by distribution function R(C, t>) which we do not specify without 
loss of information in final results. 



5 Averaged transition probabilities. 

Averaged values of objects that are local in z(t) are obtained by simple integra- 
tion with weighting function P(z,t) from (|l2|). Functional in consideration 
W nm are nonlocal in z(t), but they have a special form which allows to re- 
duce the averaging procedure to solving some parabolic differential equation^!. 
Namely, one can write down the following representation for averaged transi- 
tion probabilities at arbitrary time t 

<|c„ m | 2 )= J dCR{(,t>) J dzH nm (z)Q nm (z,t), (13) 

Functions H nm (z) are obtained from specific form of generating function for 
coefficients c mn . Functions Q n m(z,t) include ( as a parameter and are the 
solutions of the following problem 

dQ 

III 



dt 



(L - V nm )Q 



Qnm(z,t) >S(z-(), Q nm {z,t) ► 0, (14) 

*-►*> l|z|H°o 

where V nm = p nm z z , p 00 = p Q1 = 1, p 1Q = p n = 3. Operator L is defined 
in (f[2|). Formulas (|b3|)-(|l4|) are exact and give probabilities W nm at the limit 
t — ► +00. 

Let's suppose that random forces fi and f 2 act with a constant power after 
switching on and then are switched off at the moment t e . Let's also assume that 
t e is large enough to allow the replacement of Q n m(z, t e ) by the stationary limit 
Qnm(z) = um Qnm(z, t) . In such case one can obtain the following expression 

t — >+oo 

for averaged probabilities W nm 

w nm = J dbdt 2 dz 3 Q% m (t; 1 ,b,&)H„ m (a, 6, 6) , (is) 

where function Q^ m (zi, z 2 , Z3) satisfies the stationary problem 
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and there are representations for the first several functions H nm : 
#00(6,6,6) = ~ ~ ~ /^y^ exPI ^77- * [6(6 +Mi) -6 -M2] 



(16) 



Mh)\y/^j I £(6) 



#oi(6,6,6) = 



2// 



/ii = -d 5 + \l Tj— (di cos (3 + d 2 sin/?) 



2v 



s(6) 



+2 



(W) 



M2 = -de + \i~FT~ (^3 cos/? + <i 4 sin/3) 

'in 



(d 2 + d|)(l + p) - 2^/p (d 2 - dj) cos 5 + 2d x d 2 sin<5 



3 + 



(d 2 d 4 + dids)(l + p) + 1\[p [d\d\ + d 2 d%) sin 6 + (did?, ~ d 2 di) cos<5 



6+ 



+ 



(di + dj)(l + p) + 2^/p (c?4 - d|) cos 5 - 2d 3 d 4 sin<5 



(17) 



where 5 = Si + S 2 . For regular functions 6(0 an d *7o(0, defined by equations 

6 + (*)6 = o, 6 (t) ~ e in ^ , 6 (*) = 6i (t) + *6 2 (t) . 

i— » — oo 

r/ + Og(t)?7o = i-bC*), Vo(-oo) = ?)o(-oo) = 
we have used the following representations 

~ Cie 40 ""** + C 2 e- <no «**, Ci = \C\\e tSl , C 2 = \C 2 \e l8 \ 



i 



% (t) = __(^ 0< r + ^d), d(t) 

v ^ 

Also we have used designations 



\/2f2i 



^(t')F (t')dt' 



Co 



Ci 



, d = lim d(t) = V^e 4/3 . 



dl = £,Ol{t e ), d 2 = £o2(*e); d 3 = £oi(*e)> <^4 = £o2(*e), 
d5 = J?o(*e), d 6 = Vo(t e ). 

Obtained formula ( |l5| ) gives an approximate value for W nm . It is necessary 
to emphasise that it is true only for finite £i, £2, while the reducing t\ or £2 leads 
to increasing the time (t e — 1>) necessary for setting the stationary distribution. 
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